arXiv: 1504.02425vl [math.DS] 9 Apr 2015 


SHILNIKOV PROBLEM IN FILIPPOV 
DYNAMICAL SYSTEMS 

DOUGLAS D. NOVAES^ AND MARCO A. TEIXEIRA^ 


Abstract. In this paper we introduce the concept of a sliding Shilnikov 
orbit for 3D Filippov systems. Versions of the Shilnikov’s Theorems are 
provided for those systems. Specifically, we show that arbitrarily close 
to a sliding Shilnikov orbit there exist countable infinitely many sliding 
periodic orbits, and for a particular system having this kind of connection 
we investigate the existence of continuous systems close to it having an 
ordinary Shilnikov homoclinic orbit. We also prove that, in general, a sliding 
Shilnikov orbit is a co-dimension 1 phenomena. Furthermore we provide a 
family of piecewise linear vector fields as a prototype of systems having 
a sliding Shilnikov orbit. 


1. Introduction and statement oe the main results 

The study of nonsmooth dynamical systems produces interesting and amaz¬ 
ing mathematical challenges and plays an important part of so many appli¬ 
cations in several branches of science (see, for instance, na nm [B E] and 
the references therein). The present work focuses on the analysis of a typical 
phenomenon that occurs in this area which evidences a striking resemblance 
to Shilnikov homoclnic loop 

Consider a smooth three dimensional vector held for which p G is a hy¬ 
perbolic saddle-focus equilibrium admitting a two dimensional stable manifold 
and an one dimensional unstable manifold. In the classical theory of dynam¬ 
ical systems a Shilnikov homoclinic orbit T of this vector held is a trajectory 
connecting p to itself, bi-asymptotically. Under suitable genericity conditions 
this connection is a codimension one scenario, and its unfolding depends on 
the saddle quantity cr = A“ -b Re(Aj 2 ), where A“ > 0, Aj 2 ^ C (Re(Aj 2 ) < 0) 
are the eigenvalues of p. In this case, when a > 0 an intricate behaviour occurs 
when T is broken. Indeed, it is proved that there exists a compact hyperbolic 
invariant chaotic set S which contains countable inhnitely many periodic orbits 
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of saddle type in any sufficiently small neighbourhood of r(see, for instance, 

[m [la [ig [m uni [m u). 

In the theory of nonsmooth dynamical systems the notion of solutions of a 
discontinuous differential equation is stated by the Filippov’s convention (see 
[5]). In this context there exist some special points that must be distinguished 
and treated as typical singularities, one of those is a pseudo-equilibrium, which 
we shall introduce it formally later on this paper. This kind of singularity gives 
rise to the dehnition of the sliding homoclinic orbit, that is a trajectory, in the 
Filippov sense, connecting a pseudo-equilibrium to itself in an inhnity time at 
least by one side (future or past). Particularly a sliding Shilnikov orbit is a 
sliding homoclinic orbit connecting a hyperbolic pseudo saddle-focus pseudo 
to it self. 


The main goal of this paper is to produce versions of the Shilnikov’s The¬ 
orems for systems having a sliding Shilnikov orbit. More especially, we show, 
in Subsection 1^, that arbitrarily close to a sliding Shilnikov orbit there exist 
countable inhnitely many sliding periodic orbits, and for a particular system 
having this kind of connection we investigate, in Section the existence of 
continuous systems close to it having an ordinary Shilnikov homoclinic orbit. 
In Subsection |1.2[ we also prove that, in general, a sliding Shilnikov orbit is 


a co-dimension 1 phenomena. Furthermore, in Section we provide a family 
of discontinuous piecewise linear vector helds as a prototype of systems 
having a sliding Shilnikov orbit. 


1.1. Setting the problem. In this subsection the basic theory of nonsmooth 
dynamical systems is given in order to dehne the sliding Shilnikov orbits and 
to state our main results. 

Let U be an open bounded subset of We denote by C^(f/,M^) the set 
of all vector fields X : U ^ endowed with the topology induced by 
the norm ||X||r = sup{||ZlW(a;)|| : a: G 17, i G {0, 1, .. . ,r}}. Here is the 
identity operator for r = 0, and the rth-derivative for r > 0. In order to keep 
the uniqueness property of the trajectories of vector fields in C^{U, M^) we shall 
assume, additionally, that these vector helds are Lipschitz. 

Given h : 77 —)■ M a differentiable function having 0 as a regular value we 
denote by Gj((77, M^) the space of piecewise vector helds Z in M” such that 

h{x) > 0, 
h{x) > 0, 

with X, F G C''(77, M^). As usual, system ([^ is denoted by F = (X, F) and 
the switching surface h“^(0) by S . So we are taking Gj)(X, M^) = C^{K, M^) x 


( 1 ) 


Z{x) = 


X(a;), if 
F(a:), if 
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M^) endowed with the product topology. When the context is clear we 
shall refer the sets M^) and only by and C'', respectively. 

The points on E where both vectors helds X and Y simultaneously point 
outward or inward from E dehne, respectively, the escaping E® or sliding E® 
regions, and the interior of its complement in E dehnes the crossing region E'^. 
The complementary of the union of those regions constitute by the tangency 
points between X oi Y with E (see Figure [^. 



Figure 1. Dehnition of the vector held on E following Filip¬ 
pov’s convention in the sewing, escaping, and sliding regions, 
respectively. This hgure has been gotten from [2] 


The points in E'’ satisfy Xh{^) -Yhl^) > 0, where Xh denote the derivative 
of the function h in the direction of the vector X, i.e. Xh{^) = {Vh{^), X{^)). 
The points in E* (resp. E®) satisfy Xh{^) < 0 and Yh{^) > 0 (resp. Xh{^) > 0 
and Yh{^) < 0). Finally, the tangency points of X (resp. Y) satisfy Xh{^) = 0 
(resp. Yh{^) = 0). 

Now we dehne the sliding vector field 


( 2 ) 


ZiO 


YhiQXiO - XhiQYiO 
Yh{0-Xh{0 


The local trajectory of the discontinuous piecewise diherential system x = 
Z{x) passing through a point p G is given by the Filippov convention (see 

BE!): 

(z) for p G such that h{p) > 0 (resp. h{p) < 0) and taking the origin 
of time at p, the trajectory is dehned as (pz{t,p) = Pxit,p) (resp. 
Vz{t,p) = <py(t,p)) for t G Ip. 

{ii) for p E such that {Xh) {p), {Yh) (p) > 0 and taking the origin of 
time at p, the trajectory is dehned as (pzifi^p) = ipY(t,p) for t E Ip H 
{t < 0} and (pzit^p) = ipx(t,p) for t E Ip O {t > 0}. For the case 
{Xh){p), {Yh){p) < 0 the dehnition is the same reversing time; 
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{Hi) for p G E* and taking the origin of time at p, the trajectory is dehned 
as ipz{t,p) = for t G /pfljt > 0} and ipz{t,p) is either ipx{t,p) 

or ipY{t,p) or ip^{t,p) for t E Ip H {t < 0}. For the case p G E® the 
definition is the same reversing time; 

{iv) For p G SE'^ U SE'^ U SE'^ such that the dehnitions of trajectories for 
points in E in both sides of p can be extended to p and coincide, the 
orbit through p is this limiting orbit. We will call these points regular 
tangency points. 

(v) for any other point {singular tangency points) (pz{t, p) = P for all f G M; 
Here ipw denotes the flow of a vector held W. 

Remark 1. A tangency point ^ G E is called a visible fold of X {resp. Y) if 
{X)‘^h{^) > 0 {resp. {Y)‘^h{^) < 0). Analogously, reversing the inegualities, we 
define a invisible fold. Suppose thatp is a visible fold of X such thatYh{p) > 0, 
then p is an example of a regular tangency point. In this case, taking the origin 
of time at p, the trajectory passing through p is defined as (pz{t,p) = (pi{t,p) 
for t E Ip n {t < 0} and (pz{t,p) = P 2 {t,p) for t E Ip D {t > 0}, where each 
ipi, ip 2 is either px or ipy or ip^. 

A pseudo-eguilibrium is a critical point E E'^’® of the sliding vector held, 
i.e. Z{ff) = 0. When is a hyperbolic critical point of Z, it is called a 
hyperbolic pseudo-eguilibrium. Particularly if G E^ (resp. G E®) is an 

unstable (resp. stable) hyperbolic focus of Z then we call a hyperbolic 
saddle-focus pseudo-eguilibrium or just hyperbolic pesudo saddle-focus. 

In order to study the orbits of the sliding vector held it is convenient to 
dehne the (C'’) normalized sliding vector held 

(3) Z(e) = {Yh{0 - Xh{0)Z{0 = Yh{0X{0 - Xh{0Y{0, 

which has the same phase portrait of Z reversing the direction of the how in 
the escaping region. Indeed, system ([^ is obtained from (|^ through a time 
rescaling multiplying ([^ by the function Yh{^) —Xh{^) which is positive (resp. 
negative) for G E^ (resp. E E®). 

Definition 1. Let Z = {X,Y) be a piecewise continuous vector field having a 
hyperbolic pseudo saddle-focus p E {resp. p E E®). We assume that there 
exists a tangential point q E ciE® {resp. q E dY^) which is a visible fold point 
of the vector field X such that 

{j) the orbit passing through q following the sliding vector field Z converges 
to p backward in time {resp. forward in time); 

{jj) the orbit starting at q and following the vector field X spends a time 
to > 0 {resp. to < 0) to reach p. 
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So through p and q a sliding loop T is easily characterized. We call T a sliding 
Shilnikov orbit {see Figures^^for a = 0, and[^. 

Remark 2. Given Z = {X,Y) G it is worth to say that if p & is a 

fold-regular point of Z, that isp is a fold of X {resp. ofY) such thatYh{p) ^ 0 
{resp. Xh{p) ^ 0), then the sliding vector field Z (|^ is transverse to dH atp. 
A proof of this fact can be found in [15] 

1.2. Main results. In the theory of ordinary differential equations a Shilnikov 
homoclinic orbit of a 3D vector field is a co-dimension 1 phenomenon in W. 
Our first main result shows that the sliding Shilnikov is also a co-dimension 1 
phenomenon in 



Figure 2. Unfolding Z^ = {X^, Ya) of a sliding Shilnikov orbit 

rinZo = (Xo,Uo)eDF 

Theorem A. Assume that Zq = (Xo,Yo) ^ {with r > 1) has a sliding 
Shilnikov orbit Fq and let W G be a small neighbourhood of Zq. Then there 
exists a function : lU —)■ M having 0 as a regular value such that Z gW 
has a sliding Shilnikov orbit F if and only if g{Z) = 0. 

Proof. For simplicity we assume that h{x,y,z) = z, that is S = {z = 0}. Let 
Zq = (Aq, Yq) G D'” having a sliding Silnikov orbit Fq. We assume that Fq is a 
sliding loop through a pseudo-equilibrium of a focus-saddle type po G and 
a tangential point go which is a visible fold point for the vector held Xq. The 
case when po ^ would follows similarly. 

Let 7 o = Bfiqo) fl Here Ba{qo) C S is the planar ball with center at go 
and radius r. Of course 70 is a branch of the fold line contained in the boundary 
of the sliding region We remark that in the sliding region the orbit of the 
sliding vector held is always transversal to the fold line. In addition, the orbits 
of Zq through the points of 70 converge to pq in backward time. The forward 
saturation of 70 through the how of Xq meets S in a curve po in a hnite time. 
Moreover Pq G Pq. 

Let lU be a small neighborhood of Zq G So associated to each Z G W 
we can dehne similar objects: Pzyjz and pz- Clearly Z will have a sliding 
Shilnikov orbit if and only A Pz ^ Pz 
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We may assume that, in suitable local coordinate system (x,?/) around po = 
(0, 0) G fiQ is the graph of a function y = r(x). So for Z G W, pz is given 
hy y = kz{x) = oq + aix + 02 {x) with oq, oi small parameters. 

Let pz = (x*, y*) and dehne (7 : W —)■ M by g{Z) = kz{x*) — y*. Of course 
g is a. function and g{ZQ) = 0. We prove now that 0 is a regular value of g, 
that is g'{Zo) ^ 0 . 

First of all we note that g{Z) = 0 if and only if pz G pz (sliding Shilnikov 
orbit). Since 


g'^Zo) = lim 


9{Zy) -g{Zo) 


>0 V 

for any curve Z^ G O*' converging to Zq, we can take in such a way that 
Pz = (0,0) and kz{x) = v (constant). Hence g{Z^) = v and so g'{Zo) = 1. It 
concludes the proof of this lemma. □ 


Our second main result is a version of Shilnikov’s theorem for sliding Shilnikov 
orbits. 


Theorem B. Assume that Zq = (Xq, Yq) G O'” {with r > 0) has a sliding 
Shilnikov orbit Fq and let Z^ = {Xa,Ya) G be an unfolding of Zq with 
respect to Fq. Then the following statements hold: 

(a) for a = 0 every neighbourhood G C o/Fq contains countable infin¬ 
itely many sliding periodic orbits of Zq; 

{b) for every |q!| 7 ^ 0 sufficiently small there exists a neighbourhood Ga C 
o/Fo containing a finite number N{a) > 0 of sliding periodic orbits 
of Zq. Moreover N{a) —)■ 00 when a —)■ 0; 

(c) for every neighbourhood G C o/Fq there exists Ictol 7 ^ 0 sufficiently 
small such that G contains a finite number Nciao) > 0 of sliding peri¬ 
odic orbits of Za- Moreover Ncia) —)■ 00 when a —)■ 0. 

Proof. We assume that Fq is a loop through po G and go ^ The case 
Po G E® and go G 9E® would follow analogously. 

To prove statement (a) let 7 ^ = Br{qo) O cIE* and let Sr be the backward 
saturation of 7 ^ through the flow of the sliding vector held Z. The forward 
saturation of 7 ^ through the how of X meets E in a curve pr in a hnite time. 
So 

00 

Sr Fl Pr Lj, 

i=l 

where fl Jj = 0 if i 7 ^ j. The sequence of sets (/i)^i can be taken such that 
h -t {po} (see Figure |^. 

For each i = 1, 2,..., we dehne Ji as the intersection between the backward 
saturation of G through the how of X with the curve 7 ^. Clearly Jj fl = 0 
if i 7 ^ j and Ji {go}. 
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Figure 3. A schematic representation of a Shilnikov sliding 
orbit r. 



Figure 4. Periodic orbits close to a Shilnikov sliding orbit F. 


For G and ^ G let and be the flows of the sliding 

vector held Z and X, respectively. 

In what follows we dehne the applications 'ijji : 'jr ^ Ji- For ^ G 7^ there 
exists < 0 such that ^ h] and there exists tf{^) < 0 
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such that (0)6(0) e Ji. So we take 0(0 = /^(^f (0)6(0)- Note 
that O is a composition of function, being then itself a C'^ function. 

It is easy to see that each hxed point of ipi corresponds to a sliding periodic 
orbit of Z (see Figure]^. Since for each i = 1, 2,..O is a contraction. So we 
obtain a sequence (g*)^! such that g* G Ji and 'ipi{qi) = Qi- Hence we conclude 
that there exists a sequence of sliding periodic orbits of Z passing through qi. 
The proof follows just by observing that g* —?• go. 

In what follows we prove the statements ( 6 ) and (c). Firstly for |a| 7 ^ 0 
sufficiently small we build elements 7 ", S'" and /r" similarly to the elements 
7 r, Sr and fir-, respectively. 

Since the new pseudo-equilibrium is not in /i", the intersection S'" fl /i" 
has only a hnite number N{a) of disjoint sets li. Furthermore the number of 
disjoint sets N{a) in this intersection goes to inhnity when a goes to 0, and 
they converges to {p} when i ^ 00 . From here the proof of statement {b) 
follows analogously to the proof of statement (a). 

For a hxed neighbourhood G of Fq there exists |q!| 7 ^ 0 sufficiently small 
such that Pq G G n S, because Pa —t 0 when a —)■ 0, so that pa C G fl E. From 
here the proof of statement (c) follows analogously to the proof of statement 
( 6 ). □ 


2. A PIECEWISE LINEAR MODEL 


In this section we present a 2-parameter family of discontinuous piecewise 
linear dynamical system ^3 admitting a sliding Shilnikov orbit F^^^. 

For a > 0 and P > 0, consider the following discontinuous piecewise linear 
vector held. 

^ \ 


(4) 




ZaA^^y^^) = ' 


\y 

( 


— OL 

X —13 
8q! 


if z > 0, 


YaA^^y^^) = 


a 


\ 


\ 


3q; 

8a 


if z <0. 




The plane S = {z = 0} is a switching manifold for system Q, which can be 
decomposed as E = U U S® being 

E'^ = |(a:,p,0) : y> ^ j , E" = |(a:, p, 0) : p < 


8 a 


8 a 


and E^ = 0. 
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Thus p = (0,0,0) e S®. Moreover c = (/3, 3/9^/(8a), 0) is a cuspid-regular 
singularity for system (|^ (see Figure]^. 


Proposition 3. For every positive real numbers a and (3 the following state¬ 
ments hold: 


(а) the origin p = (0, 0, 0) is a hyperbolic pseudo saddle-focus of system 
Za ,0 (|^ in such way that its projection onto S is an unstable hyperbolic 
focus of the sliding vector field Za^p ([^ associated with (|^; 

(б) there exists a sliding Shilnikov orbit connecting p = (0,0,0) to 
itself passing through the fold-regular point q = {3(3/2, 3(3'^/ (8a)) (see 
Figure . 



Figure 5. A representation of the Shilnikov sliding orbit of sys¬ 
tem (|^ after bending the y-axis by doing the change of variables 
(x,y,z) = {x,y-x‘^,z). 


Proof. We compute the sliding vector held and the normalized sliding vector 
held of (|^ as 


(5) 


Za,p{x, y) 
Za,p{x, y) 


and 


Kjfj 


^ay 

3/32 (3 

ay , - X 3—y 

^’ 8a 8 ^ 


3a 


respectively. It is easy to see that (0, 0) G is a hyperbolic focus of 
Indeed, their eigenvalues are given by 




a .^/95a 
12(3 ■ 
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It implies that the origin is a hyperbolic pseudo saddle-focus of vector held 
(|^. Moreover, since Re(A^) > 0 then (0,0) is an unstable hyperbolic focus of 
the (normalized) sliding vector held (|^. 

After a change of variables and a time rescaling expressed by 
( 6 ) = * = 
respectively, the normalized sliding vector held becomes 
(7) Z = - 6 n - . 

We note that the time rescaling (|^ reverses the direction of the how of ([^. 
The fold line is given now, in (n, n) coordinates, hj i = {{u, 1) : u E M}. 

We claim that the orbit of system ([^ starting at the point (1,1) G £ is 
attracted to the focus equilibrium (0, 0) without touching the line i. Clearly, 
going back through the transformation (|^, this claim implies that the orbit of 
system ([^ starting at the point q = (3/3/2, 3/?^/(Sa)) G is attracted, now 
backward in time, to the focus ( 0 , 0 ) without touching the fold line 

To prove the claim we shall construct a compact region 77 in the u, n-plane 
that is positively invariant through the flow of the vector held ([^. To do that, 
let m{y) = —13/108 + 9//^/13 + 5%^/169, and take the curves 

Cl = {{u, 1 ) : m(l) <//<!}, 

C 2 = {(m, —2m + 3) : 1 < M < 3/2}, 

C 3 = {(m, -91/72) : m(-91/72) <u< 3/2}, 

C 4 = {(3/2,m) : -91/72 < n < 0}. 

We dehne 77 as being the compact region delimited by the curves Cj for i = 
1, 2,..., 5 (see Figure]^. After some standard computations we conclude that 
77 is positively invariant through the how of ([^. Furthermore, the vector held 
0 has at most one limit cycle (see Theorem A of i), which is hyperbolic. So 
from the positive invariance of 77, from the stability of the equilibrium (0, 0), 
and from the uniqueness of a possible limit cycle we conclude that, if this 
limit cycles exists, then it cannot be inside 77. Applying PoincarWBendixson 
theorem we conclude that the stable focus of ([^ attracts the orbits, forward 
in time, of all points in 77 without touching the line £. The claim follows by 
noting that (1,1) G 77. 

□ 
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Figure 6 . The dashed bold line represents the line £. The 
continuous bold line delimits the compact region TZ which is 
positively invariant through the flow of ([^. The bold trajectory 
is the orbit starting at (1,1) being attracted to the focus (0, 0). 

On the other hand the vector held is also linear. Thus its orbit starting 
at q is given by 


3/3 (3/3 - 2a t) (/3 + 2a t) (2/3 - 2a t)t^ 

+ UU ) -^) -TU- 
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So for = 3/3/(2a) > 0 we have that (p+(t+,g) = p. It implies that there 
exists a sliding Shilnikov orbit of Zi^a,i 3 connecting p to itself passing 
through q. 


3. Regularization 


Shilnikov showed that any smooth 3-dimensional vector held pos¬ 

sessing a hyperbolic saddle-focus p E with a 2-dimensional stable (un¬ 


stable) manifold and an 1-dimensional unstable (stable) manifold admits a 


chaotic behaviour always when its saddle quantity a = + Re(Aj 2 ) (resp. 

a = A'^ -|- Re(A“ 2 )) is positive (negative). Tresser extended the Shilnikov’s 


results for vector fields [inj and for Lipschitz continuous piecewise 
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vector fields UTl when the Shilnikov homoclinic orbit is transversal to the sets 
of non-differentiability. 

As an immediate consequence of the main result of this section we shall 
obtain that, for each positive real numbers a and /3, every neighbourhood 
U (Z ViP oi the piecewise linear model y, z), built in the previous section, 

contains a Lipschitz continuous vector held possessing a Shilnikov homoclinic 
orbit. Moreover this vector held presents a chaotic behaviour, and any neigh¬ 
bourhood of its Shilnikov homoclinic orbit contains inhnitely many periodic 
orbits. 


Theorem C. For each positive real numbers a and jd, and for ^ > 0 small 
enough, there exists a family of continuous piecewise smooth vector fields 
5-close to having the following properties for <5 > 0 small enough. 

(a) The origin is a hyperbolic saddle-focus singularity of admitting 
an 1-dimensional stable manifold Wl and a 2-dimensional unstable 
manifold Wf; 

{b) The stable and unstable manifolds intersect each other in a Shilnikov 
homoclinic orbit ^ which is 5-close to 

(c) The saddle guantity a of the origin is negative. So any neighbourhood of 
contains infinitely many periodic orbits for every 5 > 0 sufficiently 
small. 


Before proving Theorem 0 we describe the regularization process, which is 
the main too we shall use in its proof. Roughly speaking, a regularization of a 
discontinuous system Z = (df, Y) is a one-parameter family Z^ of continuous 
vector helds such that Z^ converges to the discontinuous system when 5 —)■ 0. 
The regularized system Z^ represents a class of continuous functions approxi¬ 
mated by Z as 5 —)■ 0. The Sotomayor-Teixeira method of regularization [T3] 
takes 


l + being 

Mh) ■= P{h/5), 


where 0 : M —)■ M is a continuous function which is for s G (—1,1) such 
that 0(s) = sign(s) for |s| > 1, and 0'(s) > 0 for s G (—1,1). We call 0 a 
monotonic transition function and Z^(x) the cf-regularization of Z. 

Proof of Theorem Let 0 : M —)■ M be the following monotonic transition 
function 

1 if M > 1, 
u if — 1 < M < 1, 

— 1 if u < —1. 


0(k) = < 
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and h{x, y, z) = z. We take + <^(0 , 0, A + Bx) where is the 

0-regularization of the vector held Q. Thus the differential system induced 
by for —h < 2 ; < 5, is given by 


( 8 ) 


a z 

{f3x — 3ay — 2(3‘^)z /3x + 3ay 

250 20 ’ 

{Aay — 3f3‘^)z y + {Ax + By)z {Ax + By) 

- Wa -+-2-+■*-2- 


We note that for z > 5 bhi ,/3 = ^a,i 3 + 5(0, 0, ^4 + -Bx), and for z < —5 
p + 5(0 , 0 , A + -Bx) , which are linear vector helds. 

In order to simplify the study, we take z = 5w. Thus system ([^ becomes 


X = —aw, 

13 X + 3ay + {/3x — 3ay — 2l3‘^)w 

20 ’ 

y (4ai/-302)w {Ax + By){w + 1) 

2 + —to— + ^- 2 -■ 

Here dot denotes derivative with respect to the variable t. We note that the 
origin (0, 0, 0) is the unique singularity of system ([^ for every 5 > 0. Moreover 
we can estimate their eigenvalues as 


( 9 ) y = 

6w = 




and X ^2 = 


“ ±,^ + o(S). 


120 120 


So we conclude that the origin is a hyperbolic saddle-focus singularity for every 
5 > 0 sufficiently small, which has an 1-dimensional stable manifold W| and a 
2-dimensional unstable manifold . It concludes the proof of statement (a). 

System (§, known as slow system, can be studied using singular perturba¬ 
tion methods. Doing 5 = 0 we obtain the reduced problem 


X = —aw. 


y = 
0 = 


/3 X + 3ay + (0x — 3ay — 2f3‘^)w 

20 

y {Aay - 30^)w 
2 ^ 8a 
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which is a differential equation defined on a manifold. This manifold is ob¬ 
tained as the graph through the last equality 

f 3/3^ Aay 

( 10 ) Mo = Ux,y,mo[x,y)) : a; e M, i/ < mo[x,y) = 3^2 _ 4 ^^ 

We note that (0, 0, 0) G Mq. 

Now performing the time rescaling t = 6 t we get the so called fast system 



( 11 ) 


x' = — 6 aw, 

.{/3 X + Say + {/3x — Say — 2/3‘^)w) 


y'= 5- 


2/3 


y {Aay-SP‘^)w {Ax + By){w + l) 
^ 2 ^ 8 a 2 


that we shall denote by Fs{x,y,w). Here the prime denotes derivative with 
respect to the variable r. We note that Mo is a manifold of critical points for 
system (0 when 5 = 0 , that is 


( 12 ) 


x' = 0 , 

y'= 0 , 

, ^ y {Aay - Sf3‘^)w 
^ 2 ^ 8 a 


System (12) is known as the layer problem. 

Using systems (11) and (12) it is straightforward to prove that the solution 
(^(r, 5) = 5), 5), ‘^ 3 ( 7 ’, 5)) of system ( [II| ) such that / 3 ( 0 , 5) = 1 and 

liiut^oo </ 2 (r , 5) = (0, 0, 0) can be estimated, for 5 > 0 small enough, as 

2 S(5H SftH 

= 8 a +0{6^), ip 2 {t, 6 ) = ^e 8 a +0{S^), and 


ip 3 {t, 6 ) = e 8 a + 


S/3^t ( 

45a 


9/33 


S(5H 

4a(8a + S(3H)e~ 8 a - 32e~ 4a 


V 


spH\ 

/ 


+ 0(5^ 


So the stable manifold >V| intersects the plane tc = 1 at the point 


PS = (p(0,5) = 


58a^ 58a 


3/32 ’ 3/3 


, 1 + 0 ( 5 ^). 
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For (x, I/, w) G Aio we compute 

/ 


DFo{x,y,w) = 


0 

0 

-3/^2 


0 
0 

\ Sai/ — 


0 

4q;i/ — 3/?^ 
8 a 




/ 


Since (4ai/ — 3/3^)/{8a) ^ 0 for all the points of A^o, h follows that the 
manifold AIq is a normally hyperbolic attracting manifold for Fq. So in any 
compact set of Afo we can apply the well known first Fenichel theorem (see, 
for instance, iniiHiin]), which ensures the existence of a normally hyperbolic 
attracting invariant manifold Al^ for <5 > 0 sufficiently small of system 0 
and which is known as slow manifold. The slow manifold Ads is 5-close 
to Alo, that is Al^ = {(x , ?/, m(x, i/, 5)) : m(x, ?/, 5) = mo(x, y) + 5mi(x, y)], 
where is dehned in (10). Moreover we can compute mi as 


mi{x,y) = 


l2afi‘^{Ax + B y) 48a^/?(3/5^ x + a/?^ y — 24a^?/^'' 
{day — 3/3^)^ 


(4a y — 3j3‘^)^ 


We claim that the slow manifold M.s contains the origin for 5 > 0 suf- 
hciently small. Indeed, suppose that (0, 0, 0) ^ so it is 5-close to 
because (0, 0, 0) G AIq- Since M.s is an attracting invariant manifold for 5 > 0 
sufficiently small, it must attract the origin which is contradiction because 
the origin is a singularity. Thus we conclude that (0, 0, 0) G for 5 > 0 
sufficiently small. From similar reasons the slow manifold also contains the 
unstable manifold Wf of the singularity (0, 0, 0) for 5 > 0 sufficiently small. 

We can easily check that the slow manifold AAs intersects the plane w = 1 
transversely along the curve (x, i{x, 5), 1), where 


e(x,s)^^ + s 

8a 


16a X 

3/32 


— Ax — 


3^ 

8a 


16a 


+ 0{6^ 


Now we consider the solution {x{t, 6),y{t, 6),w{t, 5)) of system 0 starting 
at a point of the slow manifold AAs- From its invariance property we know 
that w{f) = mo{x{t,6),y{t,6)) + Smi{x{t,S),y{t,S)) + C>(5^). Substituting 
this relation in the slow system (|^ we obtain the following planar differential 
system 


X 


da'^y 

day — 3/32 


+ C1(5), 


2.,2 


y = 


3fi^x + afi^y — 2da‘^y 
6/3^ — 8a(3y 


+ 0 { 6 ), 


which is topologically equivalent to the sliding vector held ([^ for 1 / < 3/3^/(4a) 
and 5 > 0 small enough. 
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Let qs = (3(3/2,i{3l3/2,6),ms{qs))- From the proof of Proposition ^ we 
know that, for <5 = 0, the orbit starting at go = (3/3/2, 3/3^/(8a), 1 ) = (9. 1 ) 
is attracted, backward in time, to the focus (0,0,0). So, for 5 > 0 sufficiently 
small, the orbit starting at 


Qs = 


3/3 3/32 6{64a^ -36Aa/3^ -9B ^ 

Y’ 24a/3 ’ 


+ 0(6^ 


is also attracted, backward in time, to the focus (0,0,0). 

Let q^ and be the points qs and ps in the variables {x, y, z) (that is z = 5w). 
The proof will follow by showing that for some branches As and Bs the flow 
of the linear system connects the points g^ and Ps for <3 > 0 sufficiently 

small. 

For z > 1 the vector held p is equal to the linear vector held Xa,p{x, y, z) + 
5 ( 0 , 0, A + Bx). Computing its solution'0(f, (5) = ('0i(f, 5),'i/’ 2 (t, 5), V' 3 (L <3)) 
such that "0(0, (3) = g^ we obtain that 

Y(C<3) = Y -at + 0{5^), 

(3/3 — 2at)(/3 + 2af) ^(6aa^ 

W, d) = -^-+ 


36Aa/3^ - 95 / 33 ) 


03(C^) = 


8a 

(3/3 — 2a 

12 


24a/3 


+ 0(6^ 


+ — 12 + 
12 ' 


32a t 

~T~ 


{<oAa - 3B/3)t‘^ - 25a 


Since the orbit 0(t, 0) reaches transversally the plane S = { 2 ; = 0} in a hnite 
time to = 3/3/(2a), we can prove that the orbit S), for <3 > 0 small enough, 
will also reach transversally the plane z = 6 in a hnite time ts- Moreover we 
can estimate ts = 3(3/{2a) + 8{32a — 9Af3‘^)/{3(3'^) + 0{8‘^). 

Let TT : ]R3 ]^2 TT-*- : m3 _i, pg |^pg projections onto the two hrst 

coordinates and onto the last coordinate, respectively. Dehne (F{A,B,8) = 
{'i(!{ts,8) —ps)/8. It is easy to see that, for every ^ > 0 sufficiently small, 
7 r-*-5(A, 5, 5) = 0 and 


nX{A,B,8) 




32a 3A(3 


35/3^ 

8a 


+ ow) ■ 


We note that 5(^0, 5o, (3o) = 0, for some Aq, Bq, and > 0, if and only 
if the vector held ([^ (for A = Aq, 5 = Bq, and 8 = ^o) admits an orbit 
connecting the points g^^ and Ps^, that is an sliding Shilnikov orbit. Since for 
A* = 40a/(9/3^) and B* = —32a‘^/{3(3^), we have that ti(F{A* , B* ,0) = 0 and 
det( 7 r 55 (A*, 5*, 0)) = —9/3^/8 7 ^ 0, then, using the implicit function Theo¬ 
rem, we conclude that, for 5 > 0 sufficiently small, there exist two branches As 
and Bs such that 7t(F{As, Bs, 8) = 0, and As —?• A* and Bs —>■ B* when 5 0. 

It concludes the proof of statement (6). 
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Finally, we compute the saddle quantity as a = —3/3^/(Sh a) + 17q;/(12/3) + 
0{S) which negative for h > 0 small enough. The proof of statement (c) follows 
by applying the classical results for Shilnikov homoclinic orbits [ISIE]. 


□ 


4. Conclusions and further directions 

In this paper we have given versions of the Shilnikov’s Theorems for systems 
having a sliding Shilnikov orbit, which is a type of connection (introduced in 
this work) in the context of nonsmooth dynamical systems. We have also 
built, explicitly, a family of discontinuous piecewise linear vector fields, 
having this kind of connection. Studying its regularization we provided a 
bridge between the classical theory and its analogous in nonsmooth dynamical 
theory. 

Higher dimension vector helds allows the existence of many other kinds of 
sliding homoclinic connections. So the study of homoclinic sliding orbits in 
higher dimension is a natural direction for further investigations. 

Finally, let Z he & piecewise vector held in admitting a sliding Shilnikov 
orbit F. We conjecture that, for all s > r, every neighbourhood U G oi Z 
contains a vector held PF G possessing a Shilnikov homoclinic orbit F. 
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